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ABSTRACT 

We analyze scattering of the transverse optic modes by spherically symmetric Polarized Nano Regions 
(PNR) in the paraelectric phase of relaxor ferroelectrics. Calculations have done in the frame of mean- 
field model early supposed by E.lolin & J.Toulouse but depolarization field effects, DF, were taken into 
account. Elementary excitations of the system are found to be of two types - Vortex (V) and Quasi Polar 
(QP). DF decreases temperature of the local QP phase transition. The local phase transition temperature 
is found to be lower (higher) for Vortex than for QP excitations for the case of small (large) size PNR. 
Therefore both QP and Vortex condensations are possible. Depolarization field hasn't effect on the 
Vortex scattering. For the case of QP scattering DF effect is essential especially in the case when total 
angular momentum j=l. Vortex and QP (at j >2) cross section of scattering I are nullified at the small 
value of momentum q of the incident wave. However QP cross section I(j=l) is saturated at q~0.01- 
0.001, and nullified only at the ultra small value of q corresponding macroscopic value of the soft mode 
wavelength. This scattering leads to the suppression of the TO spreading wave regime and allows 
qualitative understand some features of so-called waterfall phenomena observed early by the 
Brookhaven group in the neutron inelastic scattering in relaxors. 



1. INTRODUCTION 



Over the past several decades, serious progress has been made in understanding the physical 
characteristics of the disordered ferroelectrics so-called relaxors (see for example reviews [1], [2]). 
These technically important materials (PMN, PZN, KTN and other) are characterized by chemical and 
structural disordering. Early on, Burns and Dacol [3] introduced the concept of Polarized Nano Regions 
(PNR) - small ordered (polarized) regions embedded into the host crystal - as the main feature of 
relaxors. Since then, the validity of this concept has been confirmed by numerous optical, diffuse and 
inelastic neutron scattering studies (see for example [4, 5, 6]). 

An important observation was made by Gehring et al [6, 7] in their study of the Transverse Optic mode 
(soft-mode) dispersion in the relaxor PZN-8%PT in a temperature range where PNR's should certainly 
exist. These authors found that, while perfectly well defined TO peaks were observed in the neutron 
scattering spectra (constant Q-scan) for momenta q>q , these peaks disappeared for q<q and the value 
of l/q was comparable to the expected PNR size of a few tens of angstroms. In addition, the dispersion 
curves had the shape of a "waterfall". Different explanations of this phenomenon have been proposed 
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[7 -r 9]. Gehring et al [9] attributed their results to a sharp step-like increase in the TO damping for q<q . 
However, no firm consensus has been reached as to the definite origin of the waterfall phenomenon [8]. 
Another motivation for the present theoretical work are neutron scattering measurement results of the 
transverse acoustic mode (TA) dispersion and damping in the relaxor K 1 . x Ta x Nb0 3 (KTN) [10, 11, 12]. The 
shape of the TA peaks exhibit clear deviations from a standard Lorentzian form (Q-scans) at several TA 
momentum. These deviations appear to be more or less regular, not easily explained by experimental 
uncertainty and could be interpreted as evidence for resonant scattering of the TA by PNRs. 
It is now generally accepted that polarization correlations (PNR) appear at the so-called Burns 
temperature, T B , upon cooling. It has also been argued [13] that the correlations appearing at T=T B are 
dynamic and become static at a significantly lower temperature, T*, corresponding to a local distortion 
or local phase transition. T* is in fact the temperature at which elastic diffuse scattering appears in 
neutron spectra as well as the characteristic relaxor behavior or frequency dispersion of the dielectric 
constant. Consequently, the term polar nanodomains (PND) has been proposed to describe these 
polarization regions. So far, we do not have detail information about the shape of these regions, the 
nature of their boundaries (sharp or smooth), the value of the surface energy etc. 
We study PNR effect at the soft mode dynamic in the temperature interval T*<T<T B in preposition that 
the relaxor can be approximated as a two-component system. Our model is the same as we studied 
before ([14], [15]) - a spherically symmetric PNR embedded in an isotropic medium. However our new 
analysis is including previously omitted effect of the depolarizing fields (DF) at the transverse optic 
phonon (TO) scattering by PNR. We found that DF effect is very strong in our case and leads to the 
qualitative understanding of some early misunderstood feathers of waterfall phenomenon. Gehring et al 
[9] attributed waterfall results to a sharp step-like increase in the TO damping for the momentum q<q . 
TOW is corresponded to the angular momentum j>l for the case of isotropic medium. It is known that 
the cross section of the slow particles scattering, a, o~q 4 ' [17]. Therefore cross section, a, of the TO 
scattering by PNR should be small at the small momentum q and waterfall should be suppressed in this 
case. Calculation results [14, 15] demonstrate that TO scattering by PNR is really strongly increased at 
small q and nullified at, say, q<0.01 r.l.u. in agreement with general theory [17]. We found that 
depolarization field effect also leads to the result a~q 4 ' but only at the extremely small momentum, 
q<5E-7 r.l.u., corresponding to the size 1/q ~100 |am. It is happened due to the following. Depolarization 
field interact with optical mode and created corresponding displacements, which could be titled as 
Depolarization mode (DM). DM space dependence ("V" 1 inside PNR and ~r' i(i+11 ' 1 outside of PNR) is quite 
different from ~sin {qr-cot) corresponding TO wave. The incident TO wave excites charge fluctuations at 
the PNR surface. These charge fluctuations create DF interacted with optical mode displacements and 
leading to the increasing of scattering. We found that scattering cross section a corresponding angular 
harmonic j=l, ct qp (]=1), isn't nullified at very small momentum q , say q ~ 0.001. Value of a QP (j=l) is 
steadily increased parallel to the decreasing value of momentum q, saturated at the high level (plateau) 
and further nullified at the extremely small value of momentum, q~lE-6 1E-7 after passing of the 
extremely narrow resonance peak of scattering at q=q re s- Scattering cross section is crossing to the 
dependence a~q 4 at q< q res . For example value of q res =5E-7 r.l.u. and resonance energy width ~1E-12 
meV~lHz. Probably observation of similar resonances is out of limit today experimental technique, 
taking into account problems with resolution, sample quality and so on. Scattering at q~0. 01-^0.1 r.l.u. is 
very far from resonance and so no sensitive for the quality of sample. 

Therefore waterfall isn't demised at the practically important very small momentum q if DF isn't 
suppressed due to the presence of conductivity and so on reasons. 

In general problem of the DF is old and considered in any standard book concerning electric features of 
dielectric or ferroelectric materials [18]. However the interest to the DF effect is strongly increased 
during last decades due to the creation and research of materials contained many internal boundaries. 
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The effect of DF is dependent vs. conductivity, screening and so on parameters. An interesting analysis 
of these numerous questions and necessary references could be found in the paper [19]. DF is directed 
against the direction of the electric dipoles in the simple cases [18] and therefore leads to the increasing 
of the free energy. Recently it was shown [20 - 22] that DF contribution to the free energy can lead to 
the complicated structure of the ground state, instead of simple ferroelectric we could expect for 
example the presence of the toroidal momentum of polarization. 

However the complicated geometry of the polarization displacement dynamically appears even at the 
temperature which is above temperature of the local phase transition. We found that elementary 
transverse excitations of the system can be divided into two categories - Vortex (V) and quasi polar 
(QP 2 ). Vortex excitations correspond to closed or loop polarization lines of TO displacements and QP 
contain open TO polarization lines. Depolarization field has no effect at the Vortex scattering. 
In the following parts of paper we formulate the model, similar to previously proposed [14, 15], and 
provide its new solution for the boundary conditions taking into depolarizing modes. 
Our model - a spherically symmetric PNR embedded in an isotropic medium - is the simplest. 
We assume that all Hamiltonian (Lagrangian) parameters are the same outside (index 1) and inside 
(index 2) the PNR, with exception of the temperature of the phase transition (Tel, Tc2, Tcl< Tc2) and 
the energy gaps u) (l), o) (2) of the TOW soft-mode at zero momentum, q=0. We also assume that all 
parameters are constant. Mixing of transverse and longitudinal waves at the PNR surface is taken into 
account but we assume the same gap parameter for the LO and TO in order to avoid using higher 
derivatives in our equations. In the Lagrangian, we therefore include the interaction of the polarization 
with the electromagnetic field and take into account longitudinal Coulomb forces. However, we describe 
the system with a critical TO and non critical LO behavior because the LO does not usually soften in 
ferroelectrics and because Coulomb forces tend to raise its frequency at small momentum in ionic 
crystals, leading to a large energy gap Q p » u) (l), o) (2) and only a moderate (non-critical) dependence 
on temperature. LO waves are virtually excited (Q p » u) (l), ^o(2) ) but only in a thin layer near the PNR 
surface. Depolarizing field has no effect at the Vortex dynamic but is essential for the case of QP 
excitations. We show that a local instability ("phase transition in PNR") occurs at a temperature T L <Tc 2 
where the mode frequency goes to zero locally, u) L =0 (if Tc 2 >Tcj, Tc 2 is then the critical temperature for 
the case of a very large size PNR). Value of T L is different for the case of QP and Vortex excitation. 
Depolarization field effect essentially decreases value of T L for the case of QP excitation, QPT L >VT L for 
the case of small size PNR. Therefore, the QP condensation (local phase transition) is expected to occur 
before the Vortex condensation upon cooling. For the case of large size PNR , QPT L <VT L . We also found 
that local TO modes, for which value of frequency co L < (jo (1), can exist due to reflection by the host 
medium. In the last part of our paper, we calculate QP scattering by the PNR taking into account 
depolarizing field effects and find strong resonance scattering at small TO momentum as described 
above. All calculations are exact as we do not assume that the TO-TO interaction between PNR and host 
medium is small. A short discussion of the results and future perspective are presented in the conclusion 
section. 



2. MODEL 

2.1 MODEL DESCRIPTION. 



2 Similar photon excitations are called as magnetic and electric in electrodynamics [16]. 
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We consider the dynamics of a spherically symmetric PNR embedded in an isotropic medium [14, 15]. 
The total Lagrangian density L contains the contributions of the TO and LO waves, L , the electrostatic 
field, L em , and the interaction between electric field and polarization, L int . [23, 24]. 
The total Lagrangian A is written as: 

A=\Ld 3 r, L=L +L em +L^ t (1) 

L =K - n where K a =\l 2p Q (j% / dt) 2 and 77 = +(A g 1 2)div{$) 2 + bJ^ + 1 12p a% % 2 
kk =1 ' 2 (^/ ' dx k +d h 1 dx i > > kk =$ik - 1 ^ik^qq (2) 

Ant = -e\c?a<P (4) 

Here £ is the displacement in the optic mode, p the corresponding density, e* the effective charge. A , 
and B are parameters describing the optic mode contribution to the potential energy written in a 
similar form to that in elasticity theory [25] and cp is the electric scalar potential. A similar Lagrangian 
was proposed by Hopfield [26] for the description of polariton dynamics. Our Lagrangian is different 
from that proposed by Hopfield in three aspects. First, we do not take into account polariton effects. 
Second, we take into account the space dispersion of optical modes (see fl (1) first and second terms), 
and, third, we also take into account the softening of the optic mode in the last term of fl : 

2 2( \ 2 2( ) 

o)q (I) - al \T-Tc^), coq (2) = a2 \T-Tc2 j, Tc^ < Tcj, al < a2 (5) 

A phase transition is expected to occur inside the PNR at a higher temperature than in the host medium, 
T c i<T C 2- We should also note that Hopfield used the expression L H for the description of the interaction 

/ // =<2*<»V clr 4 CK - /.jnt -he* V a ( '^cx'P > 

between electric field and polarization, 

Both expressions i H and (4) lead to the same results for the case of infinitely large medium. We prefer 
apply expression (4) which is similar to the standard description of interaction between polarization and 
electric field and polarization [1, 24] and more suitable for the boundary conditions formulation. Given 
the above Lagrangian, the dynamical equations are written in the standard form: 

d I dt(SA I S(d% a 0) / dt)) = SA 1 8% a (r), £A / Sq>{r) = (6) 
From eq.2 we then get: 

p d 2 % a /dt 2 = Ao * grad a div{$) + 2Bo * V^(^) - 

-P0°>X - e% aV - A °ldHS)dS a -2Bo*\ l^dS k (7a) 

and from eqs.3 and 4: 

" V a( V a<P ~ 4™t a ) + i(y a <p-4>Kt a )dS a = (7b) 
The electrodynamics boundary condition for the tangential component of the electric field is given by: 
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V^j - nfhV^j ) = V^?2 - n(nV^2 ) ( 7c ) 



where dS is an element of the PNR surface and n is a unit vector normal to the surface. The terms 
containing dS (7a, b) lead to particular solutions satisfying the boundary conditions at the PNR interface. 
(7b, c) are the standard boundary conditions of electrodynamics. It is important to emphasize that all 
parameters are the same inside and outside the PNR with the exception of the gaps, co^l) and coo(2), 
determined by the transition temperatures T cl and T c2 respectively. 

We now consider separate independent solutions inside and outside the PNR in the form of plane waves 
with frequency oj and momentum K K T for the optic longitudinal and transverse modes respectively. 
From eqs.7a and 7b: 
Longitudinal Optic wave (LOW) 

2 2 2 2 * 

[p Q co - AoKj^ -4/3BoKj^ - P q coq ~^L e V = ^> 

iK L <p -4m 4 L = (8a) 

Solving this system of two coupled equations, we get: 

2 2 2 22 2*2 

co - AoK L I p Q -4/3BoK L I p Q - co Q ~ n p = with Q = 4x(e ) I p Q (8b) 

Because the ion plasma frequency, fl p ~15 - 75 meV »o) , is large, the softening of the LOW can be 
neglected in practice. 

Transverse Optic wave (TOW) 

[a? - BoK T 2 /p o -co 2 ]^ T = (9) 

We should note that combination of the longitudinal (8a,8b) and transverse (9) optic mode isn't enough 
for the description of the full solution of the dynamical eq-ns , depolarization field and corresponding 
displacement mode, DM, should be taken into account. DM space dependence 

("V" 1 inside PNR and ~f J<i+1> ~ 1 outside of PNR) is quite different from ~sin(gr-wt) corresponding plane TO 
wave. DF, DM features will be discussed more detailed in the section 2.3 of this paper. 
For convenience, we write the parameter A and B in the LOW and TOW dispersion, Eqs. 8b and 9, in 
terms of the speeds C b C : 

Ao= Po [c|-4/3cg ], Bo^pqcI , £->#V^b (10) 
We also take the density parameter, p , to be equal to 1. 

2.2 MODEL SOLUTION IN THE ANGULAR MOMENTUM REPRESENTATION. 

Because the effective potential energies inside or outside the spherical PNR are independent of the 
coordinates, the solution can be decomposed into independent radial and angular parts. The angular 
part is written in terms of spherical functions. The radial part can be written in terms of Bessel (Hankel) 
functions inside (outside) the PNR with argument r/C T and rK L for the TOW and LOW respectively (8a, 8b, 
9). We are especially interested in the dynamics of transverse excitations. A similar problem was 
considered in quantum electrodynamics [16]. Our case is similar but does not quite coincide with [16] 
because the longitudinal optic modes are important in our case while they are absent in [16], and 
because the optic mode energy spectrum contains a gap. At last we should take into account the 
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presence of the depolarization fields and corresponding optics wave displacements. Because we are 
interested in scattering by a spherical PNR and the total angular momentum is conserved, it is 
convenient to use [14, 15] an angular momentum technique in the following. For convenience, we utilize 
some results from Refs. [16], [17] and [27] concerning spherical vector functions. We describe the 
optic phonon polarization in terms of the vector spherical functions, x^a), u=0, ±1, a=l,2,3 (x, y, z): 



1 
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(11) 



The phonon wave field Y m contains angular and polarization ("spin") components in the total angular 

momentum j 

representation: 



r jlM 



(v,a) = I cj m i M ¥ i m (v)Z M (a% with m + p = M and v = k/k, 



m,ju 



v = 1 , / = j,j + 1 if j * 0; / = 1 if j = and M = 

where k is the momentum, V /m is a standard spherical function [16, 17] with angular momentum / and 
its projection m and d M i mlll is the Clebsch-Gordon coefficient [16]. We should emphasize that we are 
working here in the momentum (not-space) representation. Expression (12) can be written in the form 
of orthogonal normalized spherical vector functions (integration is done over the surface of a sphere of 
unit radius). 

W V) = * C lml M Y lm > J Y*jlM(v)Y rrM ,(vWl = 8 jyS^^ (13) 

We distinguish several solutions: one longitudinal, with spherical vector V'^jMiv) (superscript -1), 
directed along the momentum k, and two transverse, V%wM ar| d (superscript and +1), with 

different parity). Several useful expressions concerning spherical wave functions [16] are used in the 
paper: 



r-u (0) (i) 

y{?)(v)*x=0, Y^(v)= 



.[kV k ]7 7M (v) 



JM 
(1) 

v (v) 

(-1) 



/( / I) ' 

k 

" 77^T7 V kV (v) ' 



(14) 



JM 
(1) 



(v) 



V 2 y+l j,j-lM 
j 



(v) 



7+1 



V 2 y+l j,j+lM 

7+T 



(v), 



2j+l Y j,j-hM (V) 



(15) 



Y. . . . . (v) + 
2 j+l JJ+IM 

We note that the functions Y }M (v), Y jM (v) correspond to Vortex (V) and Quasi polar (QP) excitations 
respectively. A Vortex excitation is a manifold of closed lines of optic mode displacements (similar to the 
magnetic field lines in electrodynamics), and therefore does not possess an electric dipole moment. A 
QP TOW displacement (15) contains a spatial component corresponding to the spherical symmetric S- 
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orbital state but also spin angular momentum, giving a non-zero electric dipole moment after 
integration over angles when j=l. We remind also that only the longitudinal wave exists when the total 
angular momentum j=0. 

Next, in order to fulfill the boundary conditions at the spherical PNR surface, the phonon wave field 
must be expressed in real space. This can be done by means of an expansion of the plane wave in terms 
of spherical harmonics [16]: 

* 2 
exp(/kr) = I gi(kr)Y lm (k/ k)Y [m (n), n = r/r, n =1, 



gl {kr )=(2^) 3 1 2 i l J l+{ 1 2 {kr) I Jk~r , 

\ Y 7 7M(k / k ) exp(/kr )dQ. k = g / (kr)Y y/|fj (r/r) (16) 



Ji +1/2 is an ordinary Bessel function. As usual, Bessel functions are used in expression (16) inside PNR to 

ensure a finite TOW displacement at the PNR center, r=0; outside the PNR, we should use the first 

(second) order Hankel H !+1 /2 1>A2> function for the case of an "outward wave" ("inward wave") instead of a 

Bessel function. According to our earlier notation, the wave vector outside (inside) the PNR is labeled as 

k=Kj(l) (k=KT{2)). In the real space representation, the wave fields corresponding to the ^ 0) it Jy), ^'jmM, 

and ^jmIv) components are transformed into the expressions ^ 0> jM , F 11 '^, and &~ 1] jMl respectively: 

(0) (0) (0) (0) 

F jM = 8j ikr)Y jM („) . f. ik r)Y jM (n), (17a) 



a/7'(7+1) . - , 1V (-1). , 

2 j+l [ ~ g J +1 ( k ^ Y j M (n) " 

/j 1 ' 1 } (kr) ( n ) + /j 1 ' " 1 } (kr)Yj^ } ( n ) ( 1 7b) 

(-1) j y+l (-1) 

V = [ ^ 8 j-l (kr) + ^ 8 j + l (kr)]Y jM ( ^ + 

[Sj-l(kr)-gj + l(kr)]YfJ(n^ 

(-1,-1) (-1) (-1,1) (1) 

fj (kr)Y} M ; (n) + /) (kr) y)^ ( n ) (17c) 

where we have introduced the coefficients/ ^ f - f f lrl> j, f p f for the expansion of the wave field 
components in terms of the orthogonal basis functions T sM [n), ^jMfn), r j^n). 
It is important to note that, when going from the momentum to the real-space representation, a 
transverse vortex (0) excitations converts to a transverse vortex excitation while a transverse quasi-polar 
(+1) converts to a combination of transverse (+1) and longitudinal (-1 ) excitations (17b). Hence, the 
index ±1 in the space representation and ±1 in the momentum representation do not have the same 
physical meaning. Polarization is normally defined with respect to k, i.e. in the momentum 
representation. However, the boundary conditions are applied in real space, so that the conversion from 
momentum to real space representation is necessary. The main physical characteristic of the TOW 
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(div(£;)=0) is of course invariant with respect to the type of representation used. The functions P 0) )M , 
& +1) jM are transverse in the momentum representation (F (0) j M v=0, F (+1, jM v=0), i.e. with respect to k. In the 
space presentation however, the function P 0) jM is still transverse but this time with respect to r and the 
function tf +1> jM contains a longitudinal Y < ~ 1> j M r*0 term (17b). This term describes near field effects and is 
proportional to 1/r 2 in the wave zone (kr»l). At a large distance from the center of the PNR, 

3/2 , / sm(kr-— ) 

kr ->oc, g/ (*r) » {2k) kr 2 (18) 

We refer to F%, as the "Vortex" and to as the "quasi polar", QP, components of the TOW. In a QP 
rf^jM type excitation, both transverse and longitudinal optical mode components will be present as a 
result of applying the boundary conditions (7a), because these BCs are stated for the tangential and 
normal components of the displacements and stresses, with respect to r and not k. By contrast, the 
Vortex Ft 0> j M excitation will correspond to a strictly transverse excitation and will not be mixed with a 
LOW, since it only contains a Y° jM term, as seen in expression (17a). 

The general displacement & in the longitudinal optical wave (8b) can also be written as the gradient of 
some scalar "potential" (19). 



t T = VO , with = G . (r)Y... (3, d>) (19) 
L o o J jM 



Writing £; in terms of the potential <$> in eqs. 7a and 7b, we get: 



AO + (co - co„ - Q ) / C T <D =0 (20a) 
o p L o 

and 

2 2 2 

1 Q 7 dGj(r) j( {co -C0f.-Q.jj) 

r )-^y^G j ( r )+ j(r)=0 (20b) 

r ^or or r ^ J QjJ< J 

Function Gj(r) is proportional Bessel or Hankel functions from argument K L r. The LOW displacement (21) 
contains not only a component ccY , ' 1> JM (n) that is perpendicular to surface of the PNR but also a 
component x.V ,1 \ M (n) that is parallel to it {n=r/r, n 2 =l). 

dGj dGj (-1) ./•(./• 1) (1) 
E = —pj— nY + G .VY = —^r- Y (n) + G . Y... (n) (21) 

The near field effects described by the term Y^ 1] j M (n)G/r in eq.21 rapidly decay at a large distance r 
from the PNR center ~l/r 2 . Also, we are interested by the case when the ion plasma frequency in (8b), 
O p , is large, O p »ix). Such a LOW will only be excited within a narrow layer of thickness 6R at the PNR 
surface, 3R~l/\2nK L \ (8b). For the case of a relaxor crystal, a typical PNR radius is R~4 -r-15 l.u. and 6R~1 
l.u. if fl p =30meVand C t =100 meV*l.u. 

Let's consider depolarization field and corresponding depolarization mode displacement. We have the 
full solution of Eqn. (7b) outside of PNR surface taking into account eq-s (19, 20a, 20b) 
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A<p - Am (VJ a ) = 0, $ L = Vd> , 0> o = Gj (r)Y jM (9, <P), 

(p=4ne (®o+Sfi), AJ,p=0, 8(p=8(p j {x)Y j M {9 ,(j>), 
2 dSfij(r) 

W (r dr ) - 7(7 + (r) = (22) 

Depolarization field 6~cpj(r) which is limited inside and outside PNR is found from Eq-n (22): 

(2) j (1) -7-1 

8m . (r) = A , r - r) + i . , r 9{r - R), 9{x < 0) = 0, 9{x > 0) = 1 (23) 
7 t> t>,-l 

where >A s , A s,-i 1} -arbitrary constants and R- PNR radius. DF leads to the additional OW displacement. 

Such depolarization mode (DM) displacements are calculated by the direct substitution DF 5<p (23) to 

the eq-n (7a): 

(2) j & 2 p (2) 

% T = V (C. _,r Y ... ), C. „ = — 2 2 ^ c , r < R (24a) 

^L,a « v 2D jM'' 27) (« z -^(2)) 5 y ' 

(1) -(7+1) Q J (1) 

% = V (C, n r y... ), C, _ = — 2 o J «, . , r > R (24b) 

« v ID yM ; ' ID (« z -^(l)) S,-l' v ; 

We note that expressions (24a, 24b) are independent vs. parameters Ao, Bo (7a) defined OW space 
dispersion. DM displacements are oscillated over time but its space dependence ("V 1 inside PNR and 
~ r -i<j+i)-i j outside f pnr) i s quite different from the expressions ~sin(gr-wt), cos(qr-wt) corresponding 
plane TO wave. DM displacements create no charge density fluctuations due to the existence of the 

(2) (1) 

equalitiesy ss = v St = o- Simultaneously DM displacements are longitudinal in a 

a L,a a L,a 

(2) (1) 

meaning curl (^^ ) = curl(<%^ ) = ■ DM displacement corresponding j=l is depended vs. spherical 

angles but independent over the value of radius r inside of PNR. The electrostatic potential q> is included 
in Eqn. (7a) in the form d$/dr a and therefore doesn't affect the TOW dynamics. Due to the equality 
A6c))=0, the electrostatic depolarizing potential component, 6$, doesn't affect the LOW dynamics (8b). 
DM displacements are excited due to the scattering of the "ordinary" excitations at the boundaries. 
We now use the general solution of the model to describe the scattering of an incident transverse optic 
(TO) plane wave by a spherical PNR. Expanding the plane wave in terms of spherical harmonics [17], we 
obtain an expression for the partial wave c/ 1 ; and total, </ o; , <t w scattering cross-section: 

cj ( ' = £ (2j+l)sin(^°' 1) ) 2 2^//t 2 (25) 
7^1 J 

s(®^ and (jj 1 ) are the phase shift for scattering of vortex and quasi-polar excitations respectively. 



2.3 BOUNDARY CONDITIONS FOR THE CASE OF VORTEX EXCITATIONS 



TOW Vortex excitations create no charge density fluctuations and no mixing with LOW at the PNR 
surface. Depolarization fields are absent in this case. Therefore boundary conditions for the case of 
Vortex excitations are the same which were found early [14, 15]: 

9 



Vfo) = 0, ^=0 (26a) 

where the auxiliary function MJ(g) represents the difference between the TO amplitudes outside and 
inside the PNR: 

VQfl = A™ffun + A*lSi*«t ~ * M (2)ff/ (26b) 

A (0> (2), A l0> in , A o Ut are the Vortex amplitudes respectively inside the PNR, incident (in) and reflected 
(out) from the PNR surface on the outside; gp g jj0Ut/ g j/in - spherical Bessel function, first and second order 
spherical Hankel function, respectively. Eq-ns (26a) describe continuing of the Vortex displacement and 
it derivative over radius at the PNR surface. 

2.4 BOUNDARY CONDITIONS FOR THE CASE OF QUASI POLAR EXCITATIONS 
Electrostatic boundary condition (7c) including expression 

so that the value Gfi6(p~j (22, 23, 24a, 24b) should be continuing at the PNR surface, that is: 

(2) 2 2 2 j 

A 2L G j W + W "V 2))/ V = 
(1) 2 2 2 -7-1 

A \L°j ( R ) + C lD (°> -a) Q (l))/n R (27) 

A 2 l -LO amplitude inside PNR, A 1L - LO out wave, ~exp(l*KL*r), amplitude outside of PNR, we remind 
that corresponding momentum K L =I*£2 P /C L is the imaginary value. 

Let's consider boundary condition DS=0 from the eq-n (7b), 

DS = I (V aCp-Am ^ a )dS a (28a) 

LOW (8a) satisfies to the boundary condition DS=0. 

TOW doesn't excite electrostatic potential; TOW contribution to the value of DS at the internal part of 
the PNR surface is following: 

* (2) 

J -Am E a dS = -Am A^f. , 
u a 2 

where A 2 is the TOW amplitude and/ 2 ^ is the value of f j calculated for the region inside PNR. 
Depolarization field and corresponding displacement (24a, 24b) contribution to the value of DS at the 
internal part of the PNR surface is following: 

* 7-1 (2) * 7-1 (« 2 -»q(2)) 

Am jr {A g -C^^Anejr ( Q 2 ~^ C 2D 
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The similar calculation can be done for the respect of contribution to the value of DS at the external part 
of the PNR surface. Finally eq-n DS=0 could be written as following: 

(2) j-l {co 2 -co 2 {2)) 

- A„f. , + jr ( 7 1)C„ „ = 

2 J j,-l J y (Y p ' 2D 

(1) (1) -y-2 (co 2 -co 2 (\)) 

~ A \iJi in 1 " A \n,Ji nut 1 + ^ ~ ^ ( n 2 (28b) 

lin j, in, — I lout j ,out ,— 1 ID 

where d 1/n , j4 lout is the TOW in and out amplitude and / i; y //n ,-i (f 1] s ,om,-i) is the same as/ 1 '" 1 ^- calculated 
for the in(out) TOW and outside PNR. 

At last we have two boundary conditions corresponding continuing OW displacement perpendicular and 
parallel PNR radius and two boundary conditions corresponding continuing surface forces a aB n B (7a). 
It is convenient to express displacements as the form 

(1) (1) (-1) (-1) 
§ = G (r)Y jM + G (r)Y jM (29) 

For the case of TOW (we omit factors A 2 , Ai in ...) 



C?( 1 ) = G <^ > =[ -J— gJ+l ( kr )+-^±L^ -.j ( *r )]=/< 1 ' 1 } (kr) , 
2 7+1 2 7 + 1 J J 

< ? (-l) = C7f- 1 >= V ^ ( .^ 1) [- gy+ i(A:r)+ gy _i(A:r)]=/j 1 '- 1 >(kr) (30a) 
For the case of the LOW: 



For the case of the depolarizing displacements, we have: 

7 (1) 7-1 (-1) 7-1 

G (r < /?) = r , G = 4JU+l)r , G = jr 

Gj(r>R)=r~(j +l \ G (1) =-7 jU+l>~ J '~ 2 , =-(7+l)r~^~ 2 (30c) 

The boundary conditions (7a) is formulated as a continuity requirement on the surface force value a aS n 6 , 
that is continuing of the values of the force components £i and E_i perpendicular and parallel to the PNR 
radius respectively. 
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(1) (-1) 

a afi n p = Y jM,a L \ + Y jM,c?-\' 

i 1 =4[ 5G(1) + ffi G (-i)i G (i) ]) 

1 dr r r 

2 SG^ 1 ) VTU^T) (1) 2 (-1) 
£_! - C £ [ — g^- - G + 7 G ] + 

2 2V7U^1) (!) 4 (-1) 
C r [ r G --G ] (31) 

TOW contribution to I x and I^: 

2 5G/D VTU+i) i (1) 

1 r r 
LOW contribution to Ei and E_i : 

2 ac^ /( / i) (-i) i (i) 

Z 1L = C j[^^ + ^^ G ~r G 1. 

^ <9r r r 

2 2 /(/I) (1) 4 (-1) 

C r [ - r G - - G ] (32b) 

Depolarization mode displacement contribution to Si and S^: 

2 j-2 
r < R,Z {DM (2) = C T 2{j - 1 ),./(./ -l)r 

r<R, ^_i DM (2)=C^2j(j-l)rJ- 2 

r>R, X lDM (l)=-cj2(j+2) ./(./•!),• ' 3 , 

r > J?, £_ lzw (1) - C T 20 + Utf + 2)r (32c) 

Six eq-ns described boundary conditions allow to calculate scattering cross section, space distribution of 
TO, LO and depolarization mode displacements and so on. We will not demonstrate here these no 
complicated but cumbersome expressions but prefer directly show results of calculations. 

3. RESULTS of CALCULATIONS 



3.1. SOFT-MODE - TEMPERATURE DIAGRAM. 

We first plot in figure 1 borrowing from the ref. [15], the soft mode frequency, co(T), inside and outside 
the PNR for the particular values of the transition temperatures, T c2 and T cl , used in the calculations 
below. General conclusions about the TOW dynamics are then obtained in the following subsections. 
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Fig.l. Soft mode energy gap versus temperature 
for the case of a large radius PNR. 7d=2K, 
7c 2 =125 K, Grl=0.6; o2=l; C o =70. u) (l) 2 =al 2 (T- 
Tcx), uj (2) 2 =a2 2 (T-Tc 2 ). Numerical values of 
parameters were estimated for the particular 
physical systems considered and are being used 
in following calculations. T transp . is the 
temperature at which the soft-mode gaps inside 
and outside the PNR are the same and the soft- 
mode not scattered by the PNR. r transp =fol 2 rCj- 
a2 2 Tc 2 )/(al 2 -a2 2 ), T transp ~194.18 K. 1- triangle 
region where localized modes are possible; 2 - 
region where an incident TOW is scattered by 
PNR; 3 -condensed PNR phase. 
Fig.l couldn't be considered as general relaxor 
phase-diagram, it is rather model diagram 
corresponding to the simple, starting 
approximation. I limited himself to the 
temperature region T<T transp (for the simplicity). 

3.2. SYSTEM INSTABILITY DUE TO THE VORTEX AND OUASIPOLAR EXCITATIONS 

One PNR embedded in a very large host medium will be self transformed to another phase in process of 
local phase transition. It is assumed that temperature T>Tc 1 so thatthe host medium is in the high- 
temperature phase. A large PNR would undergo a phase transition at 7=7 c2 (Fig.l). A finite size PNR will 
transform at a lower temperature 7= T inst <Tc 2 . We calculated the value of temperature instability, T inst , 
by the following method. The instability occurs when the frequency u) of the PNR TOW mode goes to 
zero, that is we self limit by the case of "second order local phase transition". In general our calculations 
are similar to ref. [14, 15], however we took in account depolarization field effect. Therefore, we should 
find the non-zero solutions of boundary conditions corresponding to both w=0 and A Tlin =0. The 
boundary conditions are written in the form of a matrix equation: MA = 0, where the components of 
the vector d are the amplitudes A T2f A Tlout ,A L2 , A Llout , C2D, C1D and the 6x6 matrix M depends on two 
parameters -temperature, T, and PNR radius, R. Eq-n det(/W)=0 could have several solutions Tinst for 
both vortex and QP excitations, but only the first ones, corresponding to the highest temperature, are 
plotted in Fig. 2. 
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4 6 8 10 12 14 16 



R, l.u. 

Fig. 2. Temperature of the local instability for 
the case of Vortex, V T inst , and Quasi Polar, QP 
T/ n s f , excitations (depolarization field "in" and 
"off"), respectively, vs PNR radius /?. Total 
angular momentum j=l. Values of parameters: 
7c 3 =125 K, Tc 2 =2 K, al=0.6 meV/K, o2=l meV/K, 
Q p =20 mV, C o =70 meV/r.l.u., C L = 100 meV/r.l.u. 
These values of parameters will be used (as a 
rule) in the rest part of paper. Two QP 
instability mode contributions are shown for 
the depolarization field in case. 

Temperature of instability is smaller for larger 
total angular momentum j> 1 and smaller PNR 
radius R. This phenomenon is explained as in 
[15] by the presence of centrifugal barrier. 
Depolarization field effect is absent for the case 
of Vortex excitation. Depolarization field leads 
to the strong decreasing of the value of T inst for 
the case of QP excitations. It is explained by 
following. DF is directed against the direction of 
the electric dipol es in the simple cases [18] and 
therefore leads to the increasing of the free 
energy and decreasing of the value of T inst . 
Recently it was shown by the molecular 
dynamic calculations [20 - 22] that DF 
contribution to the free energy can lead to the 
complicated structure of the ground state, 
instead of simple ferroelectric we could expect 
for example the presence of the toroidal 
momentum of polarization. 
Suppose that some PNR with radius R appears 
in process of sample cooling. If DF effect is 
suppressed due to the presence of the impurity 
charges and so on reasons then QP 



condensation will appears in the process of 
further cooling [15]. Vortex condensation could 
exist at the lower temperature, VT instt , <QPT inst . 
Probably it is not contradicted to the most 
experimental data - PNR is described as 
polarized pieces of crystal. This sequence of 
events will be change for the case of non- 
screened DF: Vortex condensation will appears 
in the process of cooling (Fig. 2) if PNR radius 
R>8 l.u. QP condensation could exist at the 
lower temperature, VT instt/ > QPT inst . Authors' 
publication [20 - 22] supposed to observe 
Vortex condensation in the quantum dots. 
Unfortunately the definite experimental 
evidence concerning this phenomenon is still 
absent. 

It seems reasonably to note some indirect 
evidence of the Vortex condensation in KTN 
relaxor (KTai_ x Nb x 3 , x=0.1-K).15). Several 
experimental attempts undertaken at NIST and 
LLB by J.Toulouse et al (J. Toulouse, E. lolin, B. 
Hennion, D. Petitgrand, R. Erwin, G. Yong, 
unpublished) to observe external electric field 
effect at the neutron scattering by the 
transversal acoustic phonons in KTN weren't 
successful, electric field effect was absent. 
However Xu et al [29] observed strong effect of 
the electric field at TAW scattering at PZN- 
4.55% PT relaxor which is explained by 
reorientation ordering PNR by the electric field. 
It seems reasonably to suppose that effect of 
depolarizing field is more strongly suppressed 
(screened) in chemical and valence 
inhomogeneous PZN-4.55 PT than in valence 
and ion radius homogeneous "weak relaxor" 
KTN. The effect of the electric field will be 
absent in KTN if vortex condensation appears in 
this crystal. 

Decreasing of temperature leads to the 
appearing of the new mode instability (Fig. 2). 
Number of these modes is increased for the 
case of large PNR size. For example mode 
instability appears at 7=112.3, 100.2, 82.7, 52.1, 
and 11.9 K for the case R=15. 
A quasi-polar excitations show the existence of 
a non-zero electric dipole moment in the 
spherical layers but only for j=l. Following the 
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previous Eqs (17b) and (29, 30a, 30b and 30c), 
the average dipole moment density, P, 
associated with the QP excitation can be 
written: 




Here dO is spherical angles. The unit vector Xu 
(11) defines the direction of the dipole moment 
P and is fixed extrinsically by interaction with 
chemical short range ordering, dislocations and 
so on. We note that depolarization 
displacement contribution to the average dipole 
moment P isn't depended vs r inside PNR (30c) 
and is equal to zero outside of PNR. 



1.4 -, 




' 1 1 1 1 1 ' 1 ' 1 ' 1 1 1 1 1 1 T - 

-2 2 4 6 8 10 12 14 16 
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Fig. 3 Radial layer distribution of the QP average 
dipole moment P near the temperature of the 
corresponding local phase transition: mode 1, 
QP7" /nst =92.932K; mode 2, QP7" )nst =53.167 K 
(critical modes, frequency co=0). PNR radius 
R=7,j=l. 



Dipole moment P is continuous function (in contrast to the case of suppressed depolarizing field and 
very large LOW frequency, O p — »oo [14, 15]). It is interesting that for the case of mode 2 value of P has 
different sign at the PNR center and side (also in contrast with [14, 15] result). 



3.3 LOCALIZED OPTIC MODES 

Local optic modes (LM) are characterized by an infinite lifetime as determined by the condition lm(co)=0. 
They are in general discrete and exist in the "cage" created by the host media; the LM frequency is 
lower than the frequency of modes propagating outside the PNR, colm < ^o(l)- Therefore, LMs cannot be 
"directly" excited by a propagating incident TO wave and only shallow LMs can scatter it effectively. 
Quasi localized modes, QLM, have an energy that is greater than the soft mode energy gap of the host 
medium. They result from the effective long-time interaction between an incident TOW and the PNR. 
They are seen as higher amplitude of the propagating TO mode in the vicinity of the PNR and can lead to 
the resonant scattering. We limit ourselves here by the analysis of the LM. Calculations were done 
similar (in general) to the previously described [14, 15]. 
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Fig.4. Frequency u) of the Vortex and QP 
localized excitations, LM; j=l. Several high and 
low frequency branches are shown as a function 
of temperature. The PNR radius R =12, that is 
Vortex condensation appears at the 
temperature T=110K before QP condensation in 
process of cooling (Fig.2). Also shown is the 
o -i t t temperature dependence of the host medium 

100 ,,o ,20 130 140 150 ,60 ,70 ,80 190 soft mode gap, F, which sets an upper limit for 

the value of the LM frequency. 

Depolarization field has no effect at the Vortex local modes. However depolarizing field increases 
frequency of the QP LM and even creates new LM branch. General interpretation of these data is similar 
to the described in [14, 15]. Vortex and QP LMs could be studied by means of inelastic neutron 
scattering (in general). 
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3.4. TOW SCATTERING BY THE POLARIZED NANO REGION: STRONG EFFECT OF THE 

DEPOLARIZING MODES. 



The scattering cross-sections of transverse QP and V excitations have been calculated using eq. (25). 
Depolarizing modes aren't excited in process of Vortex scattering, corresponding V scattering cross 
section, I v , was calculated and analyzing before [14,15]. This V-scattering demonstrates the presence of 
the strong long wave resonances explained by the existence of the quasi-localized states [14, 15]. For 
the case of very small incident wave momentum q u the Vortex dynamics is similar to the scattering of a 
slow particle in quantum mechanics [17] and partial cross section of scattering, o v ~qi for q^O.Ol and 
y'=l. The similar results were found for the case of QP scattering in approximation where depolarizing 
field effect was omitted, o QP ~q* [14, 15]. 

Here we calculate TOW QP scattering taking into account depolarizing field, DF, effect (see fig. 5-10). We 
found that o QP (j=l) isn't nullified at very small momentum q t , say ql~0.001, and nullified of the 
extremely small value of momentum, ql~lE-6 + 1E-7 after passing of the extremely narrow peak of 
scattering. 




0.10 

9, 



a) 



Fig. 5a. QP scattering by PNR (total and partial 
section I). The incident TOW momentum q x is 
expressed in units n/a, parameter a is similar to 
the lattice constant; R=6, T=96K. The horizontal 
dash line marks the PNR geometrical cross- 
section. 




b) 

Fig.5b. Phase of QP scattering by PNR. R=6, 
T=96K. The resonance appears when phase = 
7i/2. Position of the cross section maximum isn't 
coincided with resonance one due to the 
presence of a factor 1/q/ in the Eq. (25). 
We should underline unusually dependence 
phase of scattering at j=l vs q lt phase ~ at 
the very small value of momentum q 3 . This 
circumstance leads to the existence of plateau 
at the cross section o(j=l) at the small 
momentum of the incident TOW (Fig.5a). 
Similar plateau is absent at the QP scattering at 
j>2 and Vortex scattering at j>l. 
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Fig.7a. QP scattering by PNR.R=6, T=150 K. 
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Fig.6a. QP scattering by PNR.R=6, T=130 K. 
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Fig.7b Phase of QP scattering by PNR. R=6, 
T=150 K. 
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Fig.6b. Phase of QP scattering by PNR. R=6, 
T=130 K. 



8a. QP scattering by PNR. R=6, T=170 K, log 
scale. 
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Fig.Sb Phase of QP scattering by PNR. R=6, 
T=170 K. 
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Fig.9a. QP scattering by PNR (total and partial 
section I). R=6, T=180 K, log scale. 
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Fig. 9b Phase of QP scattering by PNR. R=6, 
T=180 K. 
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Fig.lOa. QP scattering by PNR.R=6, T=190 K, log 
scale. 
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Fig.lOb Phase of QP scattering by PNR. R=6, 
T=190 K. 



4. DISCUSSION 



We took into account depolarization field effects in our model mean field calculations of soft mode 
dynamics in relaxors. Depolarization field leads to the strong decreasing of the value of the temperature 
of the local phase transition, T inst , for the case of quasi polar, QP, excitations. It's happed because DF is 
directed against the direction of the electric dipoles and therefore leads to the increasing of the free 
energy value. This result is in a qualitative agreement with the molecular dynamic calculations [20 - 22]. 
DF effect is absent for the case of Vortex, V, excitations. Suppose that some PNR with radius R appears 
in process of sample cooling. If DF effect is suppressed due to the presence of the impurity charges and 
so on reasons then QP condensation will appears in the process of further cooling [15]. Vortex 
condensation could exist at the lower temperature, VT insW <QPT inst . Probably it is not contradicted to 
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the most experimental data - PNR is described as polarized pieces of crystal. This sequence of events 
will be change for the case of non-screened DF: Vortex condensation will appears in the process of 
cooling (Fig. 2) if PNR radius R>8 l.u. QP condensation could exist at the lower temperature, VT msW > 

QPT in st. 

The calculated results concerning soft-mode scattering by PNR are shown in Figs. 5 - 10 and can be 
qualitatively interpreted similar to [14, 15] (with the important and later discussed exception - QP 
scattering at j=l). For the case of very small incident momentum q v the TOW dynamics is similar to the 
scattering of a slow particle in quantum mechanics [17]. The calculated values of the QP (at j>2) and 
Vortex scattering cross-section at j>l are in agreement with the general theory of scattering [17]. 
The scattering cross section I decreases at large momentum, q^O.15. 

The value of the scattering cross section has a wide maximum in the intermediate region around 
0.05 -r 0.1. The reason for the decrease of the cross section at large momentum is as follows. Scattering 
is due to the difference in the soft mode gaps (eq.5) inside and outside the PNR, 21w = ui (l) - co (2), 
which is increased at low temperature. This potential scattering is stronger at low temperature and 
decreasing for large values of q t . The effect of scattering disappears at the temperature T transp . (see Figl) 
at which the soft-mode gaps inside and outside the PNR become the same and the soft-mode is no 
longer scattered by the PNR. 

The difference between QP and Vortex scattering is essential in the resonance region due to the 
effective long time multiply interaction between excitation and PNR. QP and V cross sections are similar 
outside of this region where the corresponding cross section can be calculated in the Born 
approximation and with a plane wave basis. QP scattering resonances are explained (similar to [14, 15]) 
by the presence of quasi localized QP states. 

Let's consider more detailed QP j=l scattering by PNR. It could be noted that value of the corresponding 
cross section a(j=l) is large and saturated at the small value of q v q^O.001 (Fig. 5-10). Is it consistent 
with law of slow particle scattering o(j=l)~q 1 " in the limit q^O ? We calculated cross section, a Tay i or , at 
the ultra small value of q x by the expansion over q t coefficients of the eq-ns described boundary 
conditions and following decision of these eq-ns. Taylor (series) expansion up to 10 th order was fulfilled 
by means of software Maple 10. We found that this procedure satisfactory approximates exact 
numerical calculation up to (72=0.01 -^0.1 (Fig. 11a) and applied this approach for calculation at the 
extremely small value of q ± . 
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7=100 K, R=6): exact numerical and low ql 
taylor (series) expansion (up to q/° order) of 
the functions including to the boundary 
conditions. 
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Whole picture of scattering at the very small 
value of is shown at the Fig. lib, 11c. 
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Fig. 11a. Comparison of the values of the 
calculated QP scattering cross-section (y=l, 
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Fig. lib. 

Fig. lib. Temperature dependence of 
the whole picture of the low q! QP 
scattering, 7=100, 130, and 150 K; j=l, 
R=6, fi p =20 meV. Scattering cross 
section a is shown in the log scale. 



T,K 


Max 
phase/ft 


q! value 
corresp. 
phase 
maximum 


96 


0.5 


0.125 


130 


0.5 


0.057 


150 


0.5 


0.039 


1/U 


U.D 


u.uiy 


180 


0.0236 


0.023 


190 


0.056 


0.037 



1E7 
1000000- 
100000 
10000 
1000 



tjQPT=150 
tjQPT=170 
oQP T=180 
tjQP T=190 





Iog10(q,) 



Fig. lie Temperature dependence of the 
whole picture of the low qi QP 
scattering, 7=150, 170, 180, and 190 K; 
y=l, R=6, Qp=20 meV. 

Plateau cross section is small at T=190K, 
that is near T tra ns P ~ 194. 18K, 
approached maximum at T=170K and 
decreasing at T=100K. The reason of 
that maximum existence is following. 
The resonance QP scattering (J=l) 
appears at the small q 1 that leads to the 
stronger scattering at T=170K (see Table 
and Fig-s 5b-10b). 
Table 



We notice the presence of ultra narrow 
extremely strong peaks at the Fig lib, 11c. 
These peaks describe ultra long TO wave 
scattering by PNR. An example of similar peak is 
shown in more details at the Fig. 12. 



essential peculiarity 




q1, r.l.u. 



Fig. 12. Ultra narrow extremely strong peak 
described very long QP TO wave scattering by 
PNR: j=l, R=6, T=100 K. 
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This peak shape is approximated by the 
expression (34) contained essential peculiarity, 



CT,. 




Peak has very wide winds, it shape isn't similar 
to the well known Lorentzian or Gaussian. 
Extremely small value of the peak position 
center q lc is corresponding to the TOW wave 
length equal to the value of several mm. Cross 
section of scattering is decreased below of this 
resonance in an agreement with general theory 
[17]. 

Ultra narrow resonance peak center displaces 
to the large value of q x with increasing 
temperature. This peak is diminished at T=190K. 
Probably this extremely narrow resonance peak 
is difficult to observe, it could be destroyed by 
sample imperfections. However here we are 
interesting mainly by effect of the very far 
wings, corresponding to the value of the cross 
section to the several orders smaller than one in 
the resonance center. Therefore we could 
expect that sample imperfections effect will be 
no strong in this case. 

TOW plateau scattering parameters are 
dependent of course vs. PNR size (Fig. 13). 

The existence of the strong TOW scattering (plateau) at q^O .01 could be explained as following. The 
components of the displacement outside of PNR wfpl, wfrl, perpendicular and parallel to PNR radius 
respectively could be written in the form: 

f _ 47TiAine'>~ u} i R> ((.q 1 R) :t -q 1 R-l') ^7TiAoute' A " l i R> ((.q 1 R)' !: +q 1 R-l') C ljD V2 



Fig. 13. Temperature dependence of TOW 
plateau scattering by PNR: ^=0.001, /?=4, 6, and 
10; /=1, £2 P =20 meV. 

Effective TOW attraction to PNR is more for the 
case of large size PNR. Therefore cross section 
peak at R=6 is closer to the T transp than peak 
corresponding /?=4. Two peaks corresponding 
R=10 appearing due to the existence of 
different mode of motion. 



wfrl 



(35a) 



Here Ain, Aout, C 1D - amplitude of in, out wave and depolarization displacement respectively. We 
omitted here contribution of the LOW displacement concentrated at the narrow layer at the PNR 
surface and (for simplicity) the angular part of the expressions. 



We have by after serial expansion over q x 
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Let's suppose for a moment that we neglect depolarization displacement, that is C 1D =0. In order to have 
the finite value of the displacement at the PNR surface, we should propose that A out xA in (l+~q 1 3 ), so that 
the phase shift of the reflected wave will be ~ qi 3 , and cross section of scattering ~ qi {6 ~ 2) ~ qi in 
agreement with general theory of scattering [17]. However C 1D *0 and displacements wfpl, wfrl will be 
finite if C 1D ~1/ qi, A out =A in (l+~q 1 1 ) (this is confirmed by the numerical calculation at the qi«l), phase 
shift ~ q! and cross section of scattering is independent vs q x (plateau). Dependence C 1D ~1/ qf is in 
agreement with expression (24b) and demonstrates large susceptibility of the DM displacements for the 
respect of the electric depolarization field at the small value of q v 

Appearance of the scattering cross section plateau at j=l leads to the serious sequences. Let's suppose 

i 

that C is the PNR concentration. TOW mean free pass length X_ f _ = ; mean free pass time 

= — - group velocity V ar « -^-q, and r -^^ = 1 We should require that A mfp >2n/ qi 

in the frame of ordinary picture of spreading TOW. However this picture of the spreading TOW will be 
violated at the small value of q^ q due to the presence of plateau at the curve a(qj if PNR 
concentration C>^- . For an example, if q =0.01, o=1000, C>1.6e-06. The limit of non-overlapping PNR 

2irer 

concentration C| im =l/(4n/3*(R) 3 )=l.l e-03. Therefore we will be near the region of non-spreading soft 
mode in this case. We have TOW over damping regime for the case Tt " /p " < i. This over damping 

condition is different from non-spreading TOW condition. Therefore we can have non-over damping 
localized or quasi localized soft mode. Of course the strong soft-mode scattering leads to the some 
effect on the TOW dispersion too. We hope to consider the complicated effect of dispersion in more 
details separately. 

Gehring et al [9] attributed waterfall results to a sharp step-like increase in the TO damping for the 
momentum q<q . However here could be some trouble. 

TOW is corresponded to the angular momentum j>l for the case of isotropic medium. Therefore cross 
section, a, of the TO scattering by PNR should be small at the small momentum q and waterfall should 
be suppressed in this case. Calculation results [14, 15] demonstrate that TO scattering by PNR is really 
strongly increased at small q but nullified at, say, q<0.01 r.l.u. in agreement with general theory [17]. 
Depolarization field effect also leads to the result a~q" at j=l but only at the extremely small 
momentum, q<5E-7 r.l.u., corresponding to the size 1/q ~100 urn. It is happened due to the following. 
Depolarization field interact with optical mode and created corresponding displacements, which could 
be titled as Depolarization mode (DM). DM space dependence ("V" 1 inside PNR and ~r' ja+1> ' 1 outside of 
PNR) is quite different from ~sin(qr-wt) corresponding TO wave. The incident TO wave excites charge 
fluctuations at the PNR surface. These charge fluctuations create DF interacted with optical mode 
displacements and leading to the increasing of scattering. Scattering cross section a corresponding 
angular harmonic j=l is steadily increased parallel to the decreasing value of momentum q up to the 
very large resonance value a at q=q res and crossing to the dependence a~q 4 at q< q res . This resonance is 
extremely narrow. For example value of q res =5E-7 r.l.u. and resonance energy width ~1E-12 meV~lHz. 
Probably observation of similar resonances is out of limit today experimental technique, taking into 
account problems with resolution, sample quality and so on. Scattering at g~0.01-r0.2 r.l.u. is very far 
from resonance and so no sensitive for the quality of sample. 
Therefore plateau of scattering is existed and above mentioned trouble is over. 

Real waterfall was observed [9] at the momentum qi<<\ , qo~ 0.1. The value of q is much more than 
beginning of plateau of scattering calculated above. Therefore model parameters optimization and 
crystal anisotropy should be taken into account for the direct comparison with the existing experimental 
data. 
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It's almost impossible to observe electric DF by means of thermal neutron scattering. However 
depolarization mode displacements are excited due to the scattering of the "ordinary" thermal TO 
excitations at the PNR boundaries and probably could be observed by means of inelastic nuclear 
neutron scattering near the Bragg peaks in crystals. 
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